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Abstract 

A two-stage Stirling-type U-shape pulse tube cryocooler driven by a 10 kW-class linear compressor was 
designed, built and tested. A special feature of the cold head is the absence of a heat exchanger at the cold end 
of the first stage, since the intended application requires no cooling power at this intermediate temperature. 
Simulations where done using SAGE-software to find optimum operating conditions and cold head geometry. 
Flow-impedance matching was required to connect the compressor designed for 60 Hz operation to the 40 Hz 
cold head. A cooling power of 12.9 W at 25 K with an electrical input power of 4.6 kW has been achieved up 
to now. The lowest temperature reached is 13.7 K. 
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(j 1. Introduction 

>-> The recent development efforts for High- 
^ Temperature- Superconducting (HTS) applications 
i— i like transformers, generators or motors create a 
demand on large scale cryocoolers(author?) |[T). 

> While Gifford-McMahon (GM)- and oil-lubricated 
Stirling-coolers are presently the prevalent choice in 

l/-) HTS power applications(author?) Q3 |2l [3j , pulse 
1 J tube coolers (PTC) promise a higher reliability and 
O longer maintenance intervals. Large scale PTCs of 
the GM-type have been used for low temperature 
O cryocooling since several years now, mostly in the 

> 4 K region. While it is possible with this type of 
^ cooler to achieve large cooling powers in the 30 K 

^ region, they also suffer from relative short mainte- 
nance intervals of the compressor. Also the effi- 
ciency of the compressor/rotary-valve combination 
is quite low. Therefore, high frequency Stirling-type 
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PTCs with oil-free linear compressors became more 
attractive since they promise better cooling perfor- 
mance than GM-type coolers as well as longer main- 
tenance intervals. During the last years, several large 
scale single-stage Stirling type PTCs for operating 
temperatures between 60 and 80 K have been suc- 
cessfully developed(author?) [ffl|5][6]|. 

The work presented in the present paper is part 
of our ongoing efforts to develop a Stirling type 
pulse tube cryocooler for neon recondensation near 
25 Kfauthor?) [0[8l|9]|. Among others, this temper- 
ature is required to operate a BSCCO-based 400 kW 
HTS motor developed at Siemens for ship propul- 
sionfauthor?) [0. The low temperature allows for 
higher current densities inside the rotor coils and thus 
for a more compact design. During our high-power 
PTC development some limitations arose from a new 
kind of regenerator streaming that results from a poor 
thermal communication of the regenerator gas in 
transverse direction, as is also described in(author?) 
ifTOl . This streaming starts when the cold end tem- 
perature falls below a critical value and/or when a 
critical acoustic power (or mass flow) is exceeded. 
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Because of an excess heat transfer to the cold end, the 
streaming significantly degrates the cooling perfor- 
mance. Up to now, the critical operating conditions 
for the onset of such a streaming cannot be quanti- 
tatively predicted(author?) IfTOl . Possible solutions 
like enhancing the radial thermal conductivity of the 
regenerator matrix imply other losses caused by heat 
transfer inside the regenerator(author?) [8J. So far 
the only solution to avoid decreasing the regenera- 
tor performance is to lower the amount of acoustic 
power. 

A two-stage, gas-coupled cryocooler provides a 
way to solve these problems. First, a part of the 
acoustic power is directed into the first stage, thus 
reducing the acoustic power in the second stage. Sec- 
ond, the small acoustic power in the second stage al- 
lows a reduced regenerator diameter, which enhances 
the radial thermal communication of the working 
gas. Both measures reduce the risk of that kind of 
regenerator streaming. We therefore decided to de- 
velop a new, two-stage gas-coupled cold head, to 
minimize the possibility of regenerator streaming. 

2. Modelling of the cold head 

2.1. Design considerations 

The cold head was modelled by maintaining 
two constraints: Firstly, to reduce the possibility 
of regenerator streaming, the acoustic power in- 
side the regenerator should not exceed a critical 
value. Because theoretical predictions of the onset of 
such streaming are still under developmentfauthor?) 
IfTOll , we relied on our experiences from former de- 
signs(author?) [8] and limited the acoustic power 
W p v to 2 kW. It was clear that the target cooling 
power Qwtai of 60 W @ 25 K could not be obtained 
with such a low power, so the only solution was to 
build several identical cold heads for parallel opera- 
tion on one compressor. The necessary number n of 
parallel cold heads was estimated from 

T c 

Qtotai = nQ c = nn — W pV , (1) 

ih - 1 c 

where Q c is the cooling power of one cold head at 
temperature T c , T h the ambient temperature, and r] 
the efficiency of the cooler. From the pulse tube de- 
sign of Nguyen et al. (author?) ffTTTl . we estimated 



the target efficiency at T c = 25 K and Th = 300 K to 
7] = 10% of Carnot relative to pV-power. Then Eq.[T] 
gives a rounded-up number of n = 4 cold heads. 

Secondly, an important factor is that the required 
pV-power of 8 kW should be delivered by the avail- 
able 10 kW compressor (Qdrive model 2S297W) 
which was designed for an optimum working fre- 
quency of 60 Hz. However, consideration of the re- 
generator losses at low temperatures and high fre- 
quency revealed that the cold heads must be driven 
at a lower frequency in order to reach high efficiency 
(see next section). Another problem of using high 
frequencies is the proper flow straightening at high 
gas velocities at low gas displacement amplitudes. 
Since the pV-power scales roughly linearly with the 
working frequency, the complete cooler could not be 
tested with the present compressor. We therefore de- 
cided to test only one of the four parallel cold heads 
and to connect a more appropriate compressor at a 
later date. 

2.2. Cooler configuration 

While two-stage PTCs with thermally coupled 
stages are able to reach low temperatures(author?) 
IPT21 [T3l , there are difficulties to transfer heat over 
large cross-sections, which are required in high- 
frequency, high power PTCs(author?) [9|. Another 
problem is the large dead volume at the cold end 
of a regenerator with uniform cross-section that is 
caused by flow straighteners and the connection to 
the pulse tube, which has to be charged by the mass 
flow from the warm end. Compared to a gas-coupled 
cooler, where the cross-section is staged, this causes 
an increase in regenerator losses. While this may be 
solved by use of a conical regenerator, a gas-coupled 
two- stage cold head was favoured for practical rea- 
sons. 

Another advantage of a gas-coupled systems 
arises, when no cooling power is required at the 
first stage, as is the case with the intended HTS- 
motor cooling. In usual two-stage pulse tube sys- 
tems, a solid flow channel connects the first-stage 
pulse tube to the joint between the two regenerators, 
where some kind of flow separator divides the gas 
flow from the first-stage regenerator. Normally, the 
connecting channel is small in order to reduce void 
volume and to allow the cold gas from the first stage 



to exchange heat with the solid. From a thermody- 
namic point of view, each heat flow Q across a tem- 
perature difference AT at a mean temperature T pro- 
duces irreversible entropy: 



which leads to a decrease in exergy 3 = H — 
T,,S (author?) 02. The loss from eq. (2) can be 
avoided by eliminating such heat transfer areas and 
by thermally decoupling the first stage pulse and re- 
generator, e.g. by the use of a connecting line and 
flow straighteners with low thermal conductance. In 
such a system, the first stage will be purely "enthalpy 
coupled" avoiding irreversible entropy production by 
heat transfer. 

At the beginning of the design, a basic 
SAGE(author?) [15] model without compressor was 
set up to find the optimum operating conditions. Us- 
ing Sage, the regenerator mesh of the first stage was 
fixed to stainless steel with mesh #400 (25 um wire 
diameter) and the second stage was filled with lead 
spheres. The phase shifters consist of single iner- 
tance lines and a second-inlet at the second stage. 
The filling pressure was set to 2.5 MPa and the 
maximum pressure amplitude to 0.4 MPa as this is 
the limit of the 2S297W compressor. Simulations 
were performed to analyze the effect of operating 
frequency and input power. Fig. [TJ shows the cold 
head efficiency in a frequency range between 30 and 
60 Hz with a pV-power of 2 kW. In each simula- 
tion, the lengths and diameters of the regenerators, 
pulse tubes and phase shifters where optimized to 
give maximum cooling power at 25 K. Also the lead 
sphere diameter was optimized, while the geometry 
of the heat exchangers was not varied in order to 
speed up the optimization process. It can be seen 
from Fig{TJ that the efficiency increases when go- 
ing to lower frequencies. Most of the losses arise 
from regenerator inefficiency due to friction loss and 
small thermal penetration depth. While operation at 
a lower frequency seems to be better, one must also 
take into account that the compressor must be able to 
deliver the required acoustic power at low frequency. 
Therefore, as a compromise we have chosen a work- 
ing frequency of 40 Hz. 
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Figure 1: Calculated efficiency relative to pV-power 
and Carnot in a frequency range from 30 to 60 Hz 
with a pV-power of 2 kW. 

3. System setup 

Fig. [2] shows the system setup. The data acqui- 
sition setup is similar to that described in(author?) 
BE). The cold head is connected to the compres- 
sor via a transfer line, which consists of two tubes 
with different diameters connected in series. As de- 
scribed in section 4, the transfer line also acts as 
an impedance matching device. The aftercooler is 
of shell and tube type. The first stage regenerator 
is filled with #325/28 um and #400/25 um stainless 
steel screens. The flow distributor between the two 
regenerators is filled with an alternating stack of #80 
and #400 screens. The second stage regenerator is 
filled with lead spheres of about 100 um diameter. 
The lengths (diameters) of the first and second stage 
pulse tubes are 70 mm (26 mm), and 107 mm (16 
mm), respectively. Like the transfer line, the iner- 
tance tubes also consist of a series of two tubes with 
different diameters. 

The linear compressor is equipped with position 
sensors X\ and X 2 . Together with pressure sensor 
Pi located in the compression space, it is possible to 
measure the acoustic power delivered by the com- 
pressor in real-time. Pressure sensors are also lo- 
cated at the entrance to the aftercooler (P 2 ) and at the 
hot ends of the pulse tubes (Pj and P 4 ). Temperature 
sensors are located at the cold ends of the pulse tubes 




Figure 2: Schematic view of the system setup. X\\ X 2 : position sensors, P1-4: pressure sensors, 7V 5 ; T c \\ T c2 : 
temperature sensors, AC: aftercooler, Regl; Reg2: regenerators, Pti; Pti: pulse tubes, In\, In 2 : inertance 
and buffer, V\ : second-inlet. 



( T c i and T c2 ) and pairwise with a 180° azimuthal 
shift at the mid positions of the regenerators (T\, T 2 
and r 4 , r 5 ) to detect a possible flow inhomogeneity. 
There is also an additional sensor installed at the flow 
distributor (T 3 ). All temperature sensors except T c2 
are of PT-100 type. T c2 is a cernox sensor. An elec- 
trical heater is installed at the bottom of the second 
stage heat exchanger to enable the measurement of 
cooling power. The second stage is equipped with an 
adjustable second-inlet (V\). 

4. Properties of the impedance matching line 

After finishing the basic design of the cold head, 
the compressor model was added to the Sage model. 
The standard way to adapt the "60 Hz" compressor 
to the cold head and operate them at 40 Hz would 
have been to add a rather large void volume in be- 
tween them(author?) |[T6ll . However, it was found 
from our simulations that the required amplitude of 
the pressure wave cannot be achieved by the use of 
a void volume. Instead, we have chosen a trans- 
fer line, which acts simultaneously as an impedance 
matching device and as pressure wave amplifier. The 
length of the line is 6 m, which roughly corresponds 



to one fourth of the acoustic wave length of helium 
at 40 Hz. Fig. [3] shows a photograph of the setup. 
In order to reduce losses, the transfer line consists of 
two parts with different diameters of 43 mm and 30 
mm. Fig. [4] shows the calculated pressure and mass 
flow amplitudes in the transfer line. The input and 
output impedance are also given in the graph. It is 
seen that the transfer line amplifies the pressure am- 
plitude of 60 kPa at the inlet to 290 kPa at the outlet 
(cold head side). The amplification of the pressure 
amplitude is accompanied by a 38% loss of acoustic 
power between compressor and cold head, as seen 
from Fig. [5] During the validation of the computer 
model, significant differences were found between 
the measured and the calculated pressure phase shifts 
at both ends of the transfer line, which also affected 
the pV-power. We attribute this to the bended struc- 
ture of the transfer line, which cannot be considered 
in the model. Upon lengthening the transfer line in 
the model by 10% a better matching with the exper- 
imental results has been accomplished. The differ- 
ence in pressure amplitudes between the simulation 
and the measurement went down from 35 kPa to 10 
kPa (see also Fig. [9]), while the pressure phase angle 
error decreased from over 20° to less than 5° inside 




Figure 3: Photograph of the pulse tube cooler; the 
cold head is connected to the linear compressor via 
the 6 m transfer line. 
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Figure 4: Calculated pressure and mass flow am- 
plitude at different transfer line positions, acoustic 
power at the compressor side of the transfer line: 
W pK/ „=1.9kW. 
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Figure 5: Decrease of the acoustic power W pV in the 
transfer line. 



the compressor and from 5° to less than 2° at the end 
of the transfer line. 



5. Experimental results and discussion 

During the optimization process, several modifi- 
cations to the cold head were made. In the first runs, 
a flow maldistribution in the first stage pulse tube was 
found and fixed by adding more flow straightener 
mesh at the cold end. Later optimizations concen- 
trated on the flow distributor between the two regen- 
erators. Simultaneously with the measurements, the 
simulation model was constantly improved to reflect 
the exact construction of the cold head. All measure- 
ment were performed at 2.5 MPa filling pressure and 
40 Hz operating frequency. 

Fig. [6] shows the measured temperatures at the 
cold ends of the first and second stage pulse tubes 
and at the flow distributor between the regenerators 
as function of pV-power in inertance mode (second- 
inlet closed) compared to the corresponding simula- 
tion results. There is a good agreement for the cold 
end temperatures T c \ and T C 2, but a slight mismatch 
for the temperature T 3 , which we attribute to some 
flow inhomogeneities inside the flow distributor. 

Fig. [7] shows the temperature sensor readings T\, 
T2 and T4, T 5 (see Fig. [2]) at the middle positions 
of the two regenerators as function of the pV-power. 
There is a slight difference between T\ and T2 at 
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Figure 6: Measured temperatures T c \, T c2 , and Tj, in 
inertance mode (symbols) versus pV-power and com- 
parison to simulation results (solid lines) . 
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Figure 7: Measured temperatures in the middle of 
the regenerators versus pV-power. 
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Figure 8: The dashed line shows a result from sim- 
ulation with W p vjn =3.9 kW. The solid lines to the 
experimental data are guide for the eyes. 

higher pV-powers in the first regenerator, while in 
the smaller second regenerator a constant offset ex- 
ists, which is attributed to an inexact positioning of 
the sensors. 

Fig. [8] shows the cooling powers for various pV- 
powers in second-inlet mode. Before these mea- 
surements the second-inlet was adjusted for mini- 
mum no-load temperature. After the initial measure- 
ments took place, a modified cold head was manu- 
factured, which had less dead volumes than the orig- 
inal one. As seen from Fig. [8j this improved ver- 
sion is able to produce cooling power of 12.9 W at 
25 K with 3.9 kW of pV-power, corresponding to 
2.4 kW acoustic power at the inlet to the cold head, 
because of the loss in the transfer line (see Fig. [5]). 
The electric input power to the compressor is then 
4.7 kW at 90% of the full stroke. The lowest temper- 
ature of 13.7 K is reached in the same measurement 
with a slope of the load line of 1.15 W/K. Given the 
loss in the transfer line, this corresponds to an effi- 
ciency of 5.6% relative to the acoustic power at the 
cold head and Carnot, which is only half of the tar- 

The dashed 



get efficiency estimated in section 2.1 



line shows the simulated result with 3.9 kW of pV- 
power, which gives a cooling power of 19 W at 25 
K at a much larger slope of 1 .75 W/K and a slightly 
higher no-load temperature. While the slope can be 
attributed to the smaller pressure drop in the model 
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Table 1: Calculated enthalpy flows in the sec- 
ond stage regenerator and pulse tube and net cool- 
ing power with second-inlet closed at a pV-power 
Wpvjn =3.9 kW and T c i = 25 K. The values in 
parenthesis show the simulation results with addi- 
tional pressure drop. 
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Figure 9: Pressure amplitudes in inertance mode 
at different pV-powers from measurement (symbols) 
and from simulation (solid lines). 

compared to the experiment, a higher no load tem- 
perature points to some inaccuracies of the model 
parameters which are not yet identified. 

Fig. [9] compares the measured and simulated 
pressure amplitudes at different pV-powers in iner- 
tance mode (second-inlet closed). While for the pres- 
sure amplitudes in the compressor and at the hot end 
of the regenerator the agreement between measure- 
ment and simulation is good, the pressure amplitudes 
in the pulse tubes are lower in the experiment than 
in the simulation. In case of the second pulse tube 
the difference amounts for 10%. This is significant 
as the enthalpy flow in the pulse tube roughly scales 
with the square of the pressure amplitude. There- 
fore, compared to the experiment the calculated en- 
thalpy flow in the second pulse tube should be larger 
by 20%. Table[T]shows the calculated enthalpy flows 
including conduction losses in the second stage re- 
generator and pulse tube with second-inlet closed. 
A modified model which has the additional pressure 
drop simulated by adding some local loss in the sec- 
ond stage cold heat exchanger was made. These re- 
sults are shown in parentheses in table [T] Due to the 
additional pressure drop the enthalpy flow in the re- 
generator decreases by about 5 W (21%) and the en- 
thalpy flow in the pulse tube decreases by about 1 1 W 
(33%). The difference in enthalpy flows between the 
pulse tube and the second regenerator gives the net 



Pulse tube 2 31.6(20.8) 6.4(0.8) 



cooling power Q c , which is about 6 W larger with- 
out the additional pressure drop. The reason for the 
larger measured pressure drop is not yet identified, 
but is hoped to be fixed in the next cold head version. 

6. Conclusions 

We have developed a Stirling-type two-stage 
pulse tube cryocooler system for neon recondensa- 
tion near 25 K. In order to avoid problems with 
streaming instabilities in regenerators with large 
cross-section area, the system design is based on four 
smaller PTC cold heads to be operated in parallel on 
a single linear compressor. From experiments on one 
of the PTC cold heads in combination with a 10 kW- 
class compressor it is shown that the two-stage de- 
sign successfully suppresses the onset of internal re- 
generator streaming. It is also confirmed by the mea- 
surements that a solely "enthalpy-coupled" configu- 
ration of the stages without cold heat exchanger at 
the first stage is feasible, thereby reducing the losses 
in the first stage. 

For matching the compressor to the cold head at 
40 Hz operation, a 6 m long transfer line has been 
used that amplifies the pressure amplitude but also 
reduces the available pV-power at the cold head to 
about 62% of that at the compressor. A no-load tem- 
perature of 13.7 K and a cooling power of 12.9 W at 
25 K have been achieved at an electrical input power 
of 4.7 kW, which corresponds to 3.9 kW of acoustic 
power at the compressor and to 2.4 kW of acoustic 



power at the cold head. The corresponding efficiency 
at 25 K is 5.6% relative to Carnot and pV-power at 
the cold head. 

The calculated cooling power from simulation 
amounts to 19 W at 25 K with W pVJn = 3.9 kW as 
compared to the measured 12.9 W. The lower exper- 
imental cooling power can be attributed to a rather 
large pressure drop in the second stage, as found 
from a comparison between the measured and cal- 
culated pressure amplitudes. It is hoped that in a re- 
vised version of the cold head the pressure drop can 
be reduced and thus the cooling power increased to 
near 20 W at 25 K. 

Then, by using a larger linear compressor, de- 
signed for 40 Hz operation and capable of delivering 
about 10 kW of pV-power, four parallel cold heads 
should be able to deliver 80 W of cooling power at 
25 K for neon recondensation. 
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